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ABSTRACT 

A generic assumption of ionization models of quasar absorption systems is that radiation from local 
sources is negligible compared with the cosmological background. We test this assumption and find 
that it is unlikely to hold for absorbers as rare as H I Lyman limit systems. Assuming that the 
absorption systems are gas clouds centered on sources of radiation, we derive analytic estimates for 
the cross-section weighted moments of the flux seen by the absorbers, of the impact parameter, and of 
the luminosity of the central source. In addition, we compute the corresponding medians numerically. 
^ ' For the one class of absorbers for which the flux has been measured: damped Lya systems at z « 3, 

our prediction is in excellent agreement with the observations if we assume that the absorption arises 
in clouds centered on Lyman-break galaxies. Finally, we show that if Lyman-Break galaxies dominate 
the UV background at redshift 3, then consistency between observations of the UV background, the 
CO , UV luminosity density from galaxies, and the number density of Lyman limit systems requires escape 

fractions of order 10 percent. 

Subject headings: cosmology: theory — galaxies: formation — intergalactic medium — quasars: 
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Q\ • 1. INTRODUCTION 

The intergalactic medium (IGM) is thought to be photo-ionized by the UV radiation emitted by galaxies and quasars. 
Models of quasar absorption systems generally assume that the gas is exposed only to the mean background radiation 
that pervades the IGM. However, given that this background is produced by discrete sources, it cannot be strictly 
• uniform. 

Qh' The effect of fluctuations in the UV background on the statistics of the forest of weak HI Lya absorption lines, seen 
J*^ ^ in the spectra of distant quasars, has been the subject of a large number of papers (e.g., Zuo 1992; Fardal & ShuU 
^ 1993; Croft et al. 1999, 2002; Gnedin & Hamilton 2002; Linder et al. 2003; Meiksin & White 2003, 2004; Croft 2004; 
McDonald et al. 2005). These theoretical studies found that because most of the Lya lines arise in the low-density 
. IGM, far away from sources of ionizing radiation, global statistics of the forest are insensitive to fluctuations in the 
■ UV background. However, fluctuations in the H I ionization rate may become detectable in the Lya forest at at z > 5 
(e.g., Meiksin & White 2003, 2004) and fluctuations in the He II ionization rate have already been detected B,t z > 2.5 
from comparisons of the H I and He II Lya forests (e.g., Kriss et al. 2001; Smette et al. 2002; ShuU et al. 2004). 
; ^ In contrast to the Lya forest, many metal and high column density Lya absorbers are thought to arise in the 
d extended halos of galaxies (Bahcall & Spitzer 1969), or at least in regions close to galaxies. It is therefore not obvious 
that neglecting the contributions of local sources of ionizing radiation, which is currently common practice, is justified 
for such absorption systems. Indeed, it has been argued on the basis of metal-line column density ratios that there 
is evidence that some absorbers are exposed to a radiation field that is much softer than the general background, as 
would be expected if local stellar radiation were important (e.g., Giroux & ShuU 1997; Boksenberg, Sargent, & Ranch 
2003). 

Here we show that if we assume that a certain class of absorbers resides in the extended halos^ of a population 
of sources of UV radiation (such as galaxies), then one can estimate the typical flux from local sources to which the 
absorbers are exposed using only the luminosity density of the sources and the rate of incidence (i.e., number per unit 
redshift) of the absorbers. We do not need to know the relation between the cross-section for absorption and the 
luminosity of the source. 

We derive analytic formulas to estimate the moments of the local flux, as well as of the impact parameter and 
the luminosity of the central source, and apply them to galaxies and quasars. We find that UV radiation from local 
galaxies may well be important for absorbers rarer than Lyman limit (LL) systems (in agreement with the recent work 
of Miralda-Escude 2005) and is likely far in excess of the background for absorbers as rare as damped Lya (DLA) 
systems. We therefore conclude that the results from studies that employed ionization models of iVni > 10^ cm~ 
absorbers may need to be revised. 

This paper is organized as follows. In H2.1l we state and discuss the approximations and assumptions which we use to 
estimate the local flux. Section |2 . 21 contains the derivation of the mean flux, which we show to be close to the median 
for a wide range of models in H4.2I Section |31 presents two complementary methods to compare the local flux to the 

^ By "halo" we mean a roughly spherical region centered on a source. It does not need to be virialized, or even be gravitationally bound. 



2 



Schaye 



background flux and contains some general conclusions. Using more restrictive assumptions, we derive expressions 
for the moments of the cross-section weighted flux, impact parameter, and luminosity in H4.f I In S ij5.1l and 15.21 we 
compute the contributions of galaxies and quasars to the local H I ionization rate and we compare these with the 
background for various types of absorbers in H5.3I where we also estimate the global escape fraction for H ionizing 
photons. Finally, we summarize our main conclusions in fJHl 

2. ESTIMATING THE FLUX 

It is clearly beyond our means to compute the distribution of fluxes seen by a class of absorbers from first principles. 
Such a calculation would involve modeling radiative transfer in a characteristic volume of the universe, which would 
require specification of the full phase space distribution of elements and sources as well as of the spectra emitted by 
the sources. We would then still need to decide what we mean by the flux seen by an absorber as the absorbing gas 
clouds have a finite size. It is therefore necessary to make a number of simplifying assumptions. We will present and 
discuss these assumptions in H2.1l and estimate the mean flux in H2.2I In H4.2l we will show that the median flux, which 
we can only compute by making additional assumptions, is generally close to the mean. 

2.1. Model assumptions 

Consider a certain class of absorption systems with an observed rate of incidence dAf/dz. In practice, we will define 
a class of absorbers by specifying an ion and a minimum column density A^min- To enable us to compute the flux from 
local sources (e.g., galaxies) to which the absorbers are exposed, we will make the following simplifying assumptions 
and approximations regarding absorbers with N > A^min: 

1. All absorbers reside in spherical halos which are each centered around a single source. 

2. The flux seen by the absorber is dominated by its central source. 

3. The probability that a sightline with impact parameter b relative to a source with luminosity L intersects an 
absorber residing in the halo around the source is Pabs — fcov for b < A^min) and zero otherwise [i.e., the 
absorbing halo has a finite radius R{L, A^min)]- 

4. All of the gas contributing to the absorption resides at a distance R{L, A^min) from the central source. 

5. The central source is point-like. 

6. The product /cov/osc,jv, where fcsc,N is the fraction of the emitted radiation that is able to propagate to the 
absorber, is independent of the luminosity of the source. 

We will now discuss each of these assumptions in turn. 

The assumption that the absorbers are (or reside in) spherical halos centered on single sources will clearly break 
down for absorbers with gas densities around the cosmic mean, as these are typically distributed along sheets and 
filaments. However, it is probably a reasonable approximation for gas with densities > 10 times the cosmic mean 
and/or for heavy elements whose distributions are likely concentrated towards galaxies. 

The second assumption (radiation from sources in nearby halos is negligible) is conservative. We will show in t |5.1l 
that it is likely to be a good approximation if the sources are galaxies and if the flux from the central source exceeds 
the background. However, if the sources are strongly clustered, then radiation from nearby halos could provide a 
non-negligible contribution to the ionizing flux seen by the absorbers. 

Assumption three (Pabs = for 6 > P) is probably reasonable provided that we deflne a class by its cumulative rate 
of incidence dj\f /dz{N > N^nin)- Since dj\f /dz{N > A'min) must decrease with increasing A'min, the assumption implies 
that larger columns typically arise in sightlines with smaller impact parameters and hence that f{N, z) = d^N / dzdN is 
a decreasing function of A^. Although it is possible to fabricate density distributions that would violate this assumption, 
it cannot be violated systematically because that would mean that the cosmic density of the ion in question, which is 
proportional to the integral / Nf{N, z)dN, would diverge at large N. Moreover, the assumption that column density 
typically increases with decreasing impact parameter is supported by a large number of observational studies of the 
relation between absorption systems and galaxies (e.g., Lanzetta & Bowen 1990; Bergeron & Boisse 1991; Lanzetta ct 
al. 1995; Chen et al. 1998, 2001a, 2001b; Bowen, Pettini, & Blades 2002; Penton, Stocke, & ShuU 2002; Adelberger et 
al. 2003). 

One might naively think that the assumption that column density increases with decreasing impact parameter would 
break down in the region within which the ionizing flux is dominated by the central source (i.e., the proximity zone). 
This would be a problem, as this is one of the regimes of interest here. However, the arguments given above indicate 
that this cannot generally be the case. Before discussing our remaining assumptions, wc will show that this is not a 
paradox. 

Although the neutral hydrogen column density in a sightline passing close to a source can be small if the proximity 
zone is very large (as can be the case for bright quasars), the column density is unlikely to increase with the impact 
parameter for any given source. To see this, consider that in photo-ionization equilibrium the ionization balance 
is determined by the ratio of the ionizing flux to the gas density, which is independent of the radius if the flux is 
dominated by the central source and if the gas density proflle is isothermal, which should be close to the situation of 
interest here. In general, for an optically thin gas in photo-ionization equilibrium that is irradiated by a central source 



Local radiation and QSO absorbers 



3 



we have riHi oc n^^^l (x r^"^ if nn r"^^"*"")/^, where Fhi oc r^^ is the photo-ionization rate for neutral hydrogen. 
Now consider a spherical cloud with density profile n oc r~", where n is the number density of the species of interest. 
The column density in a sightline with impact parameter b is 

m ^ I J^^^^dlnl. (1) 

For a > 1 the dominant contribution to the column density comes from / < b and N ^ n{b)b. For example, for a 
singular isothermal sphere (a = 2) we have = Trn{b)b oc b~^, which is the scaling derived by Schaye (2001a, equations 
3 and 4) based on more general arguments. Hence, for reasonable density profiles the column density will increase 
with decreasing impact parameter. 

There are two reasons why we would in general expect this trend to be even stronger for metals than for HI. First, 
since metals are produced in stars, it is natural to assume that their abundances relative to hydrogen are decreasing 
functions of R. Second, the ionization balance of many of the metals of interest will asymptote to a constant as more 
stellar radiation is added, because (ordinary) stars essentially do not emit above 4 Rydberg (Ryd) and metals cannot 
be fully photo-ionized without these energetic photons^ 

Assumption four is also conservative: assuming that all the absorbing gas resides in a thin shell with radius R 
minimizes the flux from the central source. In fact, for a typical absorber it is probably a reasonable approximation 
to assume that all absorbing gas is at the maximum possible radius R. We already showed that for reasonable density 
profiles almost all the absorption takes place in the densest part of the gas cloud that is intersected by the sightline, 
at radii slightly greater than b. The median impact parameter of random sightlines through a halo of radius R is 
(6) = R/V^, slightly smaller than R. Hence, the column densities are typically dominated by gas at a distance r R 
from the source. We therefore expect that the assumption that all of the gas that contributes to the column density 
resides at a radius R will, for a typical absorber, result in an underestimate of the flux seen by the absorber by only a 
factor of a few. 

Note that because of assumption four, the flux seen by absorbers in some finite column density interval iVmin < N < 
-^max is identical to that seen by absorbers with N > A^min- 

Assumption five (point-like sources) will be a good approximation for sightlines with 6 ^ i?s , where Rs is the radius 
of the source. For smaller impact parameters it will, however, lead to a large overestimate of the fiux. Hence, our 
model should not be applied to absorbers with R < Rs. 

Finally, assumption six (/cov/esc,Ar is independent of L) can be justified only by our ignorance. Given that we are still 
a long way off from being able to compute this factor from first principles, it would be hard to justify any particular 
dependence on L. Assuming it is constant minimizes the number of free parameters. 

The fraction of photons that is able to propagate to the absorber, /csc.w , should not be confused with the conventional 
escape fraction /csc; which is the fraction of emitted photons that is able to escape the entire halo surrounding the 
source. They are different in general because the absorbers under consideration may in fact determine the global 
escape fraction /esc- Thus, we have < /esc < fesc,N < 1 in general with /esc.iv asymptoting to /esc for sufficiently 
small A'^niin and to unity for sufficiently large A^min (i-e., for column densities comparable to those surrounding the 
sources). 

The covering factor fcov (0 < /cov < 1) is unlikely to be much smaller than unity given that absorption line systems 
typically contain multiple components and that sightlines which pass near to galaxies nearly always show absorption 
by metals and hydrogen (e.g., Bergeron & Boisse 1991; Lanzetta et al. 1995; Adelberger et al. 2003). 

2.2. The mean flux 

We are now in a position to estimate the flux from the central source to which a typical absorber is exposed. We will 
compute the mean flux F by averaging the flux seen by an absorber over all source luminosities, weighted by the total 
cross-section for absorption. Under the assumptions discussed above, a source of luminosity L provides a (proper) 
cross-section for absorption of /covTri?^ and the flux seen by the absorber is 

Lfcsc,N /„\ 

The average number of absorption systems per unit redshift and luminosity along a random sightline is 

^{z,L,N^i^) = ^z,L)U.TrR\z,L,N^i^)j^{l + zf, (3) 

where ^ is the comoving luminosity function of the sources (i.e., the number of sources per unit comoving volume and 
per unit luminosity), c is the speed of light and H is the Huljble parameter. The mean flux from the central source to 
which a class of absorbers with rate of incidence dAf /dz{z, N > A^min) is exposed is thus 



dL 



dzdL 

dz 



1 r ^2 



dzdL 



FdL (4) 



^y^/esc,iv/cov J mz, L)dL. (5) 



^ Note that this imphes that for spectra with a sharp cutoff, such as those of galaxies, it is not true that the ionization balance depends 
only on the ratio of hydrogen ionizing photons to the gas density (i.e. , the ionization parameter) . 
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The important thing to note is that_the resuh is independent of R. Thus, we do not need to specify the function R{L) 
in order to compute the mean flux F. Equation (0) is a key result of this paper, it provides a means to estimate the 
flux seen by a class of absorbers from its observed rate of incidence and the observed luminosity function of sources. 

The mean flux is proportional to the luminosity density J L^dL. If the sources are galaxies, then this has the 
convenient consequence that F is approximately proportional to the comoving star formation density because most of 
the UV radiation is emitted by young stars. 

If the luminosity function is a Schechter (1976) function, 

^L)dL = cj^^(^y e-^/^'^, (6) 

then equation jS)) becomes 

'dj\f 



F{z, N„ 



dz 



-^-^^-^/csc,jv/cov^0*i*r(2 + a), (7) 



where r(a;) is the gamma function and k = 1/(0.4 In 10) if 0* is given per magnitude and k — 1 ii (j)^ is given per InL 
(unless specified otherwise, we will hereafter assume k = 1 and omit k from the equations). If, from some reason, the 
absorbers only reside in halos around objects with Lmin < L < Lmax, then the integral should of course be computed 
using these limits. For physically sensible luminosity functions (i.e., those that yield a finite luminosity density), 
J L^dL ^ and F will thus be insensitive to the limits as long as Lmin <C <C in 



^max . 



3. COMPARING THE LOCAL FLUX TO THE BACKGROUND 



To assess the importance of local sources, we need to compare the mean flux from the central source to the background 
radiation field. If the intensity of the background radiation is known, then we can directly compare it to our estimate 
of the mean flux [equation (|SJ)] and infer the critical rate of incidence below which local radiation dominates: 



dN \ _ dJV F 

Arit Fbg 

1 C(l + Z)2 

~F^s AH{z) 



(8) 



csc.N 



f,ov J L^z,L)dL, (9) 



where i^bg is the background flux. For sufficiently large N^^i^ (i.e., small d/V /dz) /esc, at 1 and the critical rate of 
incidence can be determined relatively accurately. But if iVmin is much smaller than the column in front of the sources, 
then /esc, AT /osc and the critical rate becomes proportional to the highly uncertain global escape fraction. 

On the other hand, if we assume that the population of sources under consideration dominates the background, that 
the proper mean free path of the photons Amfp is known, and that Amfp ^ c/ H{z) (such that cosmological redshift 
and evolution may be neglected), then the background flux is given by 

Fbg = (1 + z)'A,„fp(.2)/esc J L^z, L)dL. (10) 

Substituting this into equation lO, we see that the critical rate of incidence becomes independent of both the source 
luminosity function and the intensity of the background: 



dN" \ _ 1 /csc,jv/cov c 

/crit~^mfp(^) /c.c 4i/(z)(l + z)' 
1 fesc,Nfcov 



(11) 

(12) 



Az^mfp(,z) fesc 

where Zmfp is the redshift interval corresponding to the mean free path. 

For sufficiently large A^min, /csc,7V — 1 and the critical rate of incidence is inversely proportional to the highly 
uncertain global escape fraction f^sc- But for sufficiently small Nauni fcsc.N ~^ /esc and the critical rate of incidence 
becomes independent of the escape fraction. 

Hence, for column densities similar to those in front of typical sources, using equation will probably give a more 
accurate estimate, whereas using equation (I12|) will likely yield the most accurate estimate for lower column densities. 
In general, we can measure the global escape fraction /osc by comparing the two estimates, which comes down to 
comparing equation H10|l with the measured intensity of the background. Of course, such a measurement is subject to 
the uncertainties in the luminosity density and in Amfp- Moreover, if other sources contribute to the background, then 
the inferred escape fraction should be interpreted as an upper limit. 

Equation H12|l corresponds with equation (4) of Miralda-Escude (2005) if we neglect evolution of the rate of incidence 
on timescales < ^mfp, as we have done. However, Miralda-Escude (2005) implicitly assumed fcov = 1 and /osc, at = /osc- 
The last assumption may well break down for column densities similar to those in front of the sources, such as Damped 
Lya (DLA) Systems. 

The mean free path depends on the frequency of the photons. For energies between 1 and 4 Rydbergs (i.e., between 
the ionization potentials of H I and He II) , the mean free path depends mainly on the H I column density distribution 
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of absorbers with column densities similar to those of LL systems (A'hi > 1-6 x 10^^ cm^^). Unfortunately, this is 
exactly where the distribution is most uncertain, because the H I absorption lines are on the flat part of the curve of 
growth which makes it difficult to measure accurate column densities. Following Miralda-Escude (2003), if we assume 
that the column density distribution is a power-law dhf / dNui oc N-^^ with slope r/ = 1.5, as appears to be a reasonable 
fit at lower column densities (e.g., Kim et al. 2002), then it can be shown that A^fp = All/a/tTj where All is the mean 
spacing between LL systems. This implies that z^^^ = (cW /dz)-^^ ^/tt and hence that the critical rate of incidence is 
proportional to that of LL systems: 

dN \ ^ / c!AA\ fcscMfcov ^^^^ 

which is to be considered an order of magnitude estimate because it makes strong assumptions about the shape of the 
column density distribution. Note that this equation is less general as the previous ones because it applies only to 
photons with energy above (but near) 1 Ryd. 

If we assume that /csc,7v/cov ~ /esc, which is probably reasonable for A^hi lO^^'cm^^, then we have {dM / dz)c-c\t ~ 
{dhf /dz)-i,L, in agreement with Miralda-Escude (2005). Thus, we conclude that for hydrogen, local ionizing radiation is 
generally unimportant for Nm lO^^cm"^. For Lyman limit systems local radiation is important, provided that the 
sources associated with these systems dominate the UV background (we expect this to be the case if the background 
is dominated by stars, as opposed to quasars). The rate of incidence of DLAs is much smaller than that of Lyman 
limit systems (about a factor 10 smaller at z = 3; Storrie-Lombardi et al. 1994; Storrie-Lombardi & Wolfe 2000) which 
implies that local sources dominate the H I ionization rate in DLAs if the same population also dominates the UV 
background. Note that F may exceed i^bg by a very large factor if fcsc,N S> /esc, as may well be the case for DLAs. 

4. A POWER-LAW MODEL FOR R{L) 

Although we can compute the mean flux without specifying the function R{L), the same is not true for a number 
of other interesting quantities, such as the variance of the flux and the cross-section weighted mean luminosity and 
impact parameter. It is, however, possible to obtain analytic solutions for these and other quantities if we make the 
following additional assumptions: 

7. The luminosity function is a Schechter function. 

8. The function i? is a power-law of L. 



If 



L^^ 



R = R,[ — ] (14) 



(assumption 8), then the normalization factor i?» is determined by the observed rate of incidence: 

^ = Rl^fco.<I^A{l + zfT{l + a + 2/3), (15) 
dz H 

which follows from integrating equation and the flux (eq. 0) is given by 

/ T X 1-2/3 

f = f.(j;) . (16) 

where 

47ri?2 



' (.1') 



= (^) ^^^^i^/-c^^/™v0*i*r(l + a + 2/3). (18) 

Under these assumptions we can also solve numerically for the median luminosity and thus for the median flux and 
impact parameter. But before doing so in ti4.2l let us first consider what a reasonable value of /3 might be. 

It is difficult to predict what the dependence of R on L should be, which is why it is so convenient that F is 
independent of this unknown function. In particular, it could well be different for different elements and for different 

1/3 . . 

-^min- However, a reasonable guess might be that R scales with the virial radius, R oc rvir oc Ai^/j. , which yields 
P = 1/[3(1 - 7)] if M/L oc IVr. On the other hand, since TV ~ n{R)R (see ^T^i which becomes N - n{r^ir)rl^jR for 
a singular isothermal profile, it may be that R oc r^jj. (because n(rvir) is constant) and that = 2/[3(l — 7)] is thus a 
better guess. 

Semi-analytic models of galaxy formation (e.g., Benson et al. 2000) predict that mass-to- light ratios reach a minimum 
around L,, but that the dependence of M/L on M is weak. For example, 7 ~ 0.1, in which case we would expect 
f3 « 0.37 — 0.74, is a good fit to the predictions of Benson et al. (2000) (who modeled the i?-band at z = for a ACDM 
model and a mass range M ^ 10^^ — 10^^ -^^o)- 
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Assuming that R{L) is indeed a power-law, it is possible to measure (3 by imaging fields containing quasars. Chen 
et al. (2001a, 2001b) have done just this for the S-band (the bluer bands are the most relevant here since we are 
interested in ionizing radiation) at z < 1. They found /3 = 0.5 ± 0.1 for C IV and /3 — 0.4 ± 0.1 for A^hi > lO^'' cm~^, 
in good agreement with our naive estimates. 

Note that f3 — 0.5 is special as it gives a flux that is independent of luminosity: F — F^. Hence, if /3 is indeed 
close to 0.5, as both observations and (handwavy) theoretical arguments seem to suggest, then this has the important 
implication that absorbers belonging to a given class [defined by their rate of incidence dJ\f/dz(N > A^min)/- are nearly 
all exposed to approximately the same flux. 

4.1. Moments 

Armed with the fitting functions 10 and (|14|l . it is easy to derive a number of characteristics of the absorbers. 
Substituting these functions into equation (jsj, we obtain the following expression for the mean flux, 

F = F ^(' + ") (19) 
'^*r(l + a + 2/?)- ^'^> 

The dependence of F on /3 is only apparent, because c>c R^^ oc r(l + a + 2/3). As was shown in tl2.2l the mean flux 
is independent of R{L) (and thus (3) when expressed in physical units. 
More generally, it is not difficult to show that the m-th moment of the flux is given by, 

,™m ^„, r[l + a + m + 2/3(1 -m)] 
> - ^* r(l + a + 2/3) • 

Note that (F™) — F"'' \i (3 = 0.5. Similarly to the moments of the flux, we can derive the moments of the impact 
parameter 

and the luminosity 

r(l + a + 2^+m) 
^ ' * r(l + a + 2/3) ^ ' 

Recall that all these moments are weighted by the cross-section for absorption. Putting in numbers, for (3 = 0.5 and 
a = -1.6 (-1.2) wc obtain L = 0.9L, (1.3L,) and R « 0.5i?, (0.8i?,). FigureHshows F / F^ {left), R/R^ {middle), 
and LjL^^f {right), all as a function of /3 and for two values of a as indicated in the figure. 

4.2. Medians 

Because F is a monotonic function of L (which is a consequence of our assumptions 4 f ^2.1|) and 8), the median flux 
Fmed = F(inied), and it Can be seen from equation that the median luminosity Lmcd can be obtained by solving 
the following equation 

^Fi^dL = \\ ^R^dL. (23) 

It is easy to see from equation (I16II that for (3 = 0.5 we have Fy^cd = F = F^, because in this case F is independent 
of L. The left-hand panel of Fig. ^ shows how the median and the mean flux compare for other values of /3 and for 
two values of the slope of the faint-end of the luminosity function: a — —1.6 and a = — 1.2. The results shown in this 
figure are insensitive to Lmin and Lmax provided that Lmin ^ F* ^ Fmax- The mean flux is typically greater than the 
median, but they differ by less than a factor of two for 0.35 < jS < 0.85 if a = -1.6 (0.22 < /3 < 1.05 if a = -1.2). 
Note that the power-law model for R{L) becomes unphysical for low (3 because i?* ^ for /? — > — (1 -f a)/2. 

The median impact parameters are shown in the middle panel of Fig. ^ again for two values of a as indicated in 
the figure. As was the case for the flux, the medians are typically close to the means. 

For the mean flux to be an order of magnitude greater than the median, we would require /3 > 1.6. Such a large value 
would be in conflict with the observations of Chen et al. (2001a, 2001b) and would probably imply an extremely low 
value of 7 (7 < —1.4 for f3 — 2/[3(l — 7)]). Hence, although the median flux cannot be computed without specifying 
R{L), it is within a factor of a few of the mean for plausible models. This gives us confidence that equati_on jSJ gives 
a reasonable estimate of the typical flux seen by the absorbers. Thus, we expect that in general Fmcd ~ F ~ F*. 

5. APPLICATION TO OBSERVATIONS 

The two types of sources that we will consider are galaxies and quasars. We will compute the characteristic flux at 
1 Ryd as well as the H I ionization rate for two relatively well measured luminosity functions: Lyman-break galaxies 
(LBGs) at z = 3 (Steidel et al. 1999a) and quasars at z — 2.3 (Croom et al. 2004). Both of these luminosity functions 
were derived using samples selected in the rest-frame ultraviolet (UV). Hence, they do not include highly obscured 
sources, which may provide a substantial contribution to the total rates of star formation and accretion onto black 
holes in the universe. However, since such systems would likely also be absent from surveys of quasar absorbers, which 
typically target quasars that are exceptionally bright in the rest-frame UV, these luminosity functions are exactly what 
is needed for our purposes. 
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FlG. 1. — The flux from the central galaxy seen by the absorbers (left), the impact parameter relative to the central galaxy (middle), 
and the luminosity of the central galaxy (right) are all plotted as a function of /3, the exponent of the power-law R = Rt,(L/ Lt)l^ . In 
each panel both the (cross-section weighted) mean and the median are shown for two different slopes of the faint-end of the luminosity 
function: a = —1.6 and a = —1.2. The quantities plotted are all dimensionless because the flux, impact parameter, and luminosity are 
plotted relative to the characteristic quantities Ft, and L,. While the latter is the familiar characteristic luminosity appearing in the 
Schechter luminosity function (eq. 1^]), the parameters Ft, and i?* depend on /3, the rate of incidence of t he a bsorb ers, the parameters of 
the luminosity function, and the factors fcov and /esc jv- Their values can be computed using equations 1151 and Hl^l . If the mean flux 
were plotted in physical units, rather than relative to Ft, then it would be independent of /3. Note that the model predictions become 
unphysical at small /3 because Rt — > as /3 — > — (1 -|- ce)/2. 



5.1. Galaxies 

By combining their ground-based data with data from the Hubble Deep Field (HDF), Steidel et al. (1999a) find 
that, down to at least O.IL*, the z « 3 luminosity function is well fit by a Schechter function with parameters 
a = -1.60 ±0.13, = 24.48 ±0.15 (AB magnitude), and 0* = 4.45 x IQ-^/i^Mpc"^ [for^ (r2„,r2A) = (0.3,0.7) and 
(z) — 3.04]. We note that because they relied on HDF data for L < L^,, which covers only a relatively small field, the 
faint-end slope may be subject to systematic errors due to cosmic variance. 

7^,(AB) = 24.48 ± 0.15 at z = 3.04 yields L^^, = (5.72 ± 0.79) x lO^^ h'^ ergs-^Hz"^ /-i,^ at rest-frame 

A = Aeff/(1 + z) = 1715.35 A. Starburst99 (version 4.0, Leitherer et al. 1999) predicts L^(A1715)/L^(A912) = 3.6 for a 
continuous starburst (the result converges within 10 Myr) with a Salpeter initial mass function with upper and lower 
mass cutoffs of 100 Mg and 1 Mq respectively, and a metallicity Z = 0.2^0 [Lj,(A1715)/I/j.(A912) = 4.1 for Z = Zq]. 
Using this value we obtain _L^^*(A912) = (1.59 ± 0.22) x 10^^ h^'^ ergs^-^ Hz~^ /o^c ai715- Equation lO then gives 

F.(A912)lbg,.=3 = (4.43 ± 0.90) x 10-^° ergcm-^ s-^ Hz-^ ^^^ /esc.Ar.A9i2 /dM\ ' ^^4) 

/csc,A1715 \ J 

Using the same Starburst99 model for the spectral shape, we find that this corresponds to a H I ionization rate of 

fHi,LBG,.=3 = (9.94 ±2.03) X 10-12 s-i /cov#^ (^) \ (25) 

/csc,A1715 \az J 

where (/csc.Af) is the average fraction of hydrogen ionizing photons that escapes to the absorber, weighted by the 
cross-section for ionization. Neglecting evolution in the luminosity density, the values for z — 2.3 would be 0.90 times 
those for z = 3. 

It is interesting to estimate the typical impact parameter of absorbers around LBGs. Assuming the power-law 
radius-luminosity relation (eq. d]), equation p5(l yields, 

i?*,LBG,.=3 = 1.16xl02kpc f— j [/cov^(2/3-0.4)]-l/^ (26) 

where we assumed h = 0.7. Using j3 = 0.5 gives .R « 48 kpc [dM / dzY^'^ fcm'^ ■ To put this in context, damped Lya 
systems (DLAs, A'ni.min = 2 x 10^" cm"^) have an observed rate of incidence d/V /dz = 0.20 at z « 3 (Storrie-Lombardi 
& Wolfe 2000). Hence, DLAs are likely to have impact parameters of ~ 10 kpc (see also Fynbo et al. 1999; Haehnelt, 
Steinmetz, & Ranch 2000). Comparing this with the half-light ratios of bright (L > L*) LBGs, which are typically 
about 2 kpc (0.2-0.3 arcsec; Giavalisco Steidel, & Macchctto 1996) we see that the approximation of point-like sources 
is likely to be reasonable. However, we do expect our model to overestimate the local flux for absorbers that are much 
rarer than DLAs. 

3 Steidel et al. (1999a) quote 0, = 1.6 X IQ-^/jSMpc^^ mag ^ for (Q.m,^A) = (1.0,0.0). We changed the units from per magnitude to 
per InL by dividing by 0.41nlO and we converted to (Q.m,^A) = (0.3,0.7) using the effective volumes listed in their Table 3. 
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We can also use equation (|26|l to verify the validity of our assumption that the flux from nearby sources is negligible 
compared to that from the central source. Since galaxies with L ~ L* dominate the cosmic luminosity density, we would 
expect this assumption to be a good approximation if (lH-z)'^i?ij0* ^ 1 if the galaxies were randomly distributed. From 
equation (|26() we can see that this would imply dM jdz <C 4 x 10^/covr(2/3 — 0.4). However, bright LBGs are observed 
to cluster: their inferred comoving correlation length is — hhr^ Mpc for a two-point correlation function of the form 
^ = (r/ro)~^ * (Steidel et al. 1999b). Hence, the expected number of galaxies within a distance i?, of a i»-galaxy is 
about 2.5[i?*(l + z)/ro]^^'^ higher than for a random distribut ion a nd we require dN jdz ^ 1 x 10^/covr(2/3 — 0.4) for 
the flux from nearby galaxies to be negligible. We will see in H5. 31 that for z < 4 this condition is just satisfied if the 
flux (at around 1 Ryd) from the central galaxy exceeds the background. However, given the approximate nature of 
this estimate, it is possible that we have underestimated the mean local flux by a factor of a few for classes of absorbers 
for which the flux from the central galaxy is not large compared with the background. 

There is one class of absorbers for which the intensity of the local UV radiation has been measured: DLAs with 
detectable absorption by C II*. Wolfe et al. (2004) have measured the UV flux at 1500 A (rest-frame) for the 23 
absorbers in their sample of 45 z ~ 2-4 DLAs for which they were able to obtain a positive detection of C II* 
absorption. By equating inferred cooling rates to grain photo-electric heating rates (given the measured dust-to- 
gas ratios), they find that their C II* absorbers are typically exposed to a flux that is within a factor of a few of 
-47rl0-i^-5 « 4x10-18 erg cm-2s-iHz"\ Our Starburst99 spectral template predicts Fi.(1500)/F,,(1715) = 1.12, 
which yields F^(A1500) = 2x lO-^^ ergcm^^ g-i Hz"^ /cov/csc,w.ai50o//csc,ai715 for dM jdz = 0.1 (« 0.20-23/45, where 
0.20 is the rate of incidence for 2^3 DLAs measured by Storrie-Lombardi & Wolfe 2000). Thus, /or the one class 
of absorbers for which the local UV flux has been measured, the prediction of our model is in excellent agreement with 
the observations. Measurements of the local flux based on H2 lines confirm that the UV radiation field in DLAs far 
exceeds the background (e.g., Hirashita & Ferrara 2005, Srianand et al. 2005). 

The good agreement between our prediction for the fiux to which DLAs are exposed, which relied on the assumption 
that DLAs arise in the halos of LBGs, and the measurement of Wolfe et al. (2004) suggests that DLAs and LBGs 
may be drawn from the same population of galaxies, as has already been suggested on other grounds by Schaye 
(2001b) and M0ller et al. (2002). If this is indeed the case, then we can read of the predictions of our model for the 
impact parameters and luminosities from figure ^ For f3 ~ 0.5 and z = 3 we predict a median impact parameter 
med(i?) « 0.4i?, « 14 kpc « 1.8 arcsec) (see eq. jJEl) and a median luminosity med(L) « O.IL*. 

5.2. Quasars 

Using the double power-law fit to the luminosity function of quasars in the 2dF QSO redshift survey (Groom et al. 
2004), we find / L^{L)dL = 2.67 x lO-^^L, Mpc"^ and Mb,,, = -25.77 for {n,n,^i^,h) = (0.3,0.7,0.7) and z = 2.3. 
Assuming bj = B, a power-law spectrum F^, = Fy{vij^){v / v-iA^)~°''' , where i^LL is the frequency at the hydrogen Lyman 
limit, and a spectral index a^, — 1.8 (Telfer et al. 2002) this gives Li,^*(A912) = 5.23 x 10^° ergs-^Hz-^ /esc A1333 ^^"^ 
thus a mean flux 

i^.(A912)Qso..=2.3 = 5.05 X lO'^^ ergcm-^ s-^ Hz-^ f,J-fS^E^^ ' . (27) 

/csc,A1333 \a,Z J 

For a power-law spectrum the H I ionization rate is Fhi ~ 9.51 x 10* Fy{vu^) / [i + a^) . Hence, if there exists a 
population of absorbers which all reside in halos around QSOs, then such absorbers will on average be exposed to a 
local radiation field with an ionization rate 

f HI,QSO,.=2.3 = 1.00 X 10-12 S-1 /eov4^^ ('^) ^ ■ (28) 

/osc,A1333 \ az / 

Since both F and the intensity of the UV background (UVB) radiation are proportional to the luminosity density of 
the sources, the relative contributions of the central galaxies and quasars to F are identical to their relative contributions 
to the UVB. 

However, although every quasar probably resides in a galaxy, only a few percent of galaxies host an active quasar 
(e.g., Steidel et al. 2002). Therefore, the mean fiux from local quasars may not be very relevant to quasar absorption 
studies. Studies employing non-parametric statistics (such as the median and other percentiles) and studies using 
small samples of absorbers, are unlikely to be significantly affected by local quasar radiation, unless the absorbing 
species requires the presence of a hard radiation source for its existence. It should, however, be noted that regardless 
of how small the fraction of time that a quasar is active, local radiation from quasars would be important if they 
are dormant for periods comparable to or smaller than the timescale for the establishment of ionization equilibrium. 
Although this timescale is typically very short for hydrogen (< 10^ yr), it can be much longer for species with higher 
ionization potentials. 

In the remainder of this paper we will ignore radiation from local quasars. 

5.3. Comparison with the background 

How does the flux from local sources compare to that from the extragalactic background? In section Owe found 
that, relative to the background, the local fiux is given by 

F /dJV\ rdAf\'^ 
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where the critical rate of incidence can be computed using either equation @ or H12|) . The former method requires a 
measurement of the luminosity density, the background flux, and the escape fraction to the absorber (/esc, at), while the 
latter requires measurements of the mean free path and the ratio /esc,jv//csc, and also requires one to assume that the 
sources associated with the absorbers dominate the background. As the absorber column density tends to the column 
in front of the sources, /esc, at asymptotes to unity, while for much lower absorber column densities it asymptotes to 
the global escape fraction /esc- Hence, provided the mean free path is known and that the sources associated with 
the absorbers dominate the background, the second method is probably more reliable for all but the highest column 
densities. If the luminosity density, the background, and the mean free path are all known, then we can measure the 
global escape fraction fcsc by comparing the two methods. 

We will now apply both methods to hydrogen ionizing photons. For other ions the calculation is entirely analogous 
to that for H I, except that a model of the spectral shape of the UVB is required because the rates of photo-ionization 
of heavier elements by the background have not been measured. For ions with ionization potentials i?ion < 4 Ryd we 
expect the critical rates of incidence to be similar to those of H I, but for ions with Eion > 4 Ryd, such as O VI, the 
critical rates should be much smaller because stars (with the exception of Pop. Ill) emit very little radiation at these 
energies. Note, however, that even though local stellar radiation is unlikely to be important above 4 Ryd, ionization 
models of O VI absorbers (as well as other ions with ionization potentials above 4 Ryd) generally also use constraints 
from other species whose ionization balance is affected by the addition of radiation below 4 Ryd. 

First, we will compute the critical rate of incidence using method 1 (i.e., equation (SJ . For the intensity of the UV 
background we use rHi/(10~^^ s~^) — (0.06,0.67,1.3,0.9,1.0) at z = (0,1,2,3,4). These measurements were taken 
from-* Dave & Tripp (2001) (z = 0), Scott et al. (2002) (z = 1), and Bolton et al. (2005) (z = 2-4). For redshift 3 
we have already computed the mean ionization rate due to local sources (equation I25|l. We can scale this result to 
other redshifts if we assume that the luminosity density of ionizing radiation is proportional to the cosmological star 
formation density, for which we take /9*(z = 0) = 0.1p*(z = 1) and p^{z > 1) = p^(z = 1) (e.g.. Heavens et al. 2004 
and references therein). We then find that 



dM\ ^(5^ 9^ 6, 11, 11) 

dz /^j.jj /osc,A1715,z=3 



at z = (0, 1, 2, 3, 4). (30) 



Second, we compute the critical rate of incidence using method 2, which reduces to equation l|13|) for H I, if we neglect 
the variation of the mean free path with frequency and if we assume that the H I column density distribution is of the 
form dAf /dN cx N^^'^ around the Lyman limit. Making this assumption and using {dM / dz)^^ w (0.7, 1.0, 1.3, 2.0, 4.0) 
for z = (0, 1, 2, 3, 4) (Peroux et al. 2005), we obtain 

^) =(0.3,0.4,0.6,0.9,1.8)^^^^2^^^^^ at z = (0, 1, 2, 3, 4). (31) 

dz J Jcsc 

We stress that the ratio of /-factors could be greater than unity for column densities similar to those shielding the 
sources. 

Comparing equations H30(l and H31|l . we find 

f 

— = (0.06, 0.05, 0.09, 0.08, 0.16) at z = (0, 1, 2, 3, 4). (32) 

/csc,A1715,z=3 

These measurements of the global escape fractions result from the demand that the luminosity density due to LBGs 
gives rise to the measured UV background. Because the measurements of the luminosity density, the mean free path, 
and the intensity of the UV background are all uncertain by a factor of a few or more, these measurements should be 
taken as order of magnitude estimates only. Thus, we conclude that a global escape fraction fcsc ^ 10~^/csc.ai715.2=3 
is consistent with all the measurements. If quasars dominate the UVB then the escape fraction should be smaller. 

Steidel, Pettini, & Adelberger (2001) measured Fj.(A1500)/i^^(A900) = 4.6 ± l.O from their composite LBG 
spectrum (after correction for intervening absorption) whereas our spectral template predicts an intrinsic value of 
Li,(A1500)/Li,(A900) = 4.0, implying that fcsc ^ /esc,Ai500,z=3- Although this may be consistent with equation H32|l 
within the errors (which are very poorly constrained) there is some tension between the two. However, since the 
galaxies studied by Steidel et al. (2001) were unusually blue, it is important to keep in mind that their measurements 
may not be representative for the population as a whole. 

If (/esc,Ar) /cov ~ /esc, then local ionizing radiation becomes important for H I for LL systems. Absorbers with 
A'hi > 10^^ cm~^ have rates of incidence much lower than LL systems (e.g., Peroux et al. 2005) and are thus dominated 
by local radiation. Since estimates of the local UV flux have generally not been available, the implication is that 
published results from ionization models of such systems may be significantly in error. Note that for A^hi > lO^*' cm^^ 
the absorber columns are probably comparable to those in front of the sources in which case fcsc,N may be much 
greater than /esc, further boosting the local flux. 

Comparison of the rates of incidence of metal line absorbers with the critical rates indicates that UV radiation 
from local galaxies may not be negligible for strong metal systems, such as those containing detectable absorption by 
ions in low ionization stages. For example, using the data from Boksenberg et al. (2003), who decomposed a large 

* Where appropriate the measurements have been scaled to the currently favored cosmology {Q,f,h? = 0.0224, h = 0.71; Spergel et al. 
2003) using the relation V ot Hg/i^ (e.g., Ranch et al. 1997). 
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number of high-quality quasar absorption spectra into Voigt profiles, we find^ dM jdz — 4.5 for log A'i„in(C II) = 12 and 
dM jdz = 3.4 for log A^min(Si IV) = 13, both at z « 3. These observed rates of incidence are likely to be overestimates 
of the true rates since many of the quasars in the sample of Boksenberg et al. (2003) were originally selected to contain 
DLAs. Indeed, more than half of their C II systems with logiV > 12 and nearly half of their Si IV systems with 
log > 13 are associated with DLAs. Other frequently studied metal line systems with low rates of incidence include 
Mg II systems (dM jdz k, 1 for systems with an equivalent width Wa > 0.3 A at z '--^ 1; Steidel & Sargent 1992). 

It is difficult to predict exactly how the physical properties inferred from ionization models would change if local 
radiation had been taken into account, because the effects will depend on the ions involved as well as on the rate 
of incidence of the absorbers. However, if the column density of H I is used as a constraint, then the addition of 
local radiation will generally tend to increase the inferred density and to decrease the inferred size of the absorber. 
The effects on the abundances of heavy elements are more difficult to predict and will likely be different for different 
elements. 

6. CONCLUSIONS 

We have constructed an analytic model to estimate the characteristic flux from local sources of radiation to which 
quasar absorption systems are exposed. Since many studies have shown that fluctuations in the UV background are 
unimportant for the low-column density H I Lya forest (at least at z < 5), we focused on the rarer systems with higher 
H I columns and/or detectable absorption by heavy elements. 

The most important assumptions we made are that the absorption arises in a roughly spherical gas cloud centered 
on a source of ionizing radiation, that radiation from other nearby sources is negligible, and that the probability that 
a sightline intersects a total column density A > A^min is zero beyond some radius i? (which may depend on the 
luminosity L and A^min)- The first assumption will break down for the low column density Lya forest, which is thought 
to arise in a web of sheet-like and filamentary structures, but is likely to be reasonable for the rarer systems with higher 
H I columns and/or detectable absorption by heavy elements. The second assumption is conservative. We argued it 
is likely to be a reasonable approximation for those absorbers for which we predict the flux from the central galaxy 
to exceed the background (around 1 Ryd). The third assumption implies that higher columns arise in sightlines with 
smaller impact parameters, which we argued must typically be the case. 

We showed that the mean flux from the central galaxy is proportional to the luminosity density of the sources 
and inversely proportional to the rate of incidence dM /dz{N > N^in) of the absorbers, but that it does not depend 
on the function i?(L, Amin). Assuming a power-law dependence R cx , we derived analytical expressions for the 
cross-section weighted moments of the flux, of the impact parameter, and of the luminosity. In addition, we computed 
the corresponding medians numerically. We found that for reasonable values of /?, the distribution of fluxes is narrow 
and the mean (which is independent of /?) is of the same order of magnitude as the median. For the special case of 
/3 — 0.5 all absorbers belonging to a class with a given rate of incidence are exposed to the same flux. Interestingly, 
both observations and naive theoretical arguments suggest that (3 is close to this value. This implies that there exists 
a characteristic flux to which absorbers with A^ > N^i^ are exposed, which can be estimated using equation (|SJ). 

We applied our model to two relatively well-studied sources of ionizing radiation: Lyman-break galaxies (LBGs) at 
z « 3 and quasars at z « 2. We argued that galaxies are more relevant to quasar absorption systems because quasars 
are thought to be dormant for most of their lifetimes. 

Our predictions are in excellent agreement with the observations for the one class of absorbers for which the flux 
has been measured: DLAs at z « 3 with detectable absorption by C II*, Since our calculation assumed that the 
absorbers are centered on LBGs, this agreement suggests that DLAs and LBGs may be drawn from the same underlying 
population of galaxies. We predicted (for /3 0.5) that DLAs at z « 3 typically have impact parameters of order 
10 kpc and that the median luminosity of their LBG counterparts is about O.IL*. Predictions for classes of absorbers 
with other rates of incidence can easily be obtained from the expressions derived in §fJ21and21 

We used two different methods to compare the mean, local flux to the background. The first takes the luminosity 
density and the escape fraction to the absorber as inputs, and requires knowledge of the background intensity. The 
second takes the mean free path for photons and the ratio of the escape fraction to the absorber and the global 
escape fraction as inputs, and relies on the assumption that the sources associated with the absorbers dominate the 
background. We found that consistency between the two methods, which comes down to the requirement that the 
observed sources make up the observed background, requires a global escape fraction of order 10 percent. 

For absorbers as rare or rarer than Lyman limit systems (which have dM /dz « 2 at z = 3), local H ionizing 
radiation becomes important. (Sub-)DLA systems are likely dominated by local radiation. Studies that have modeled 
the ionization balance of such systems assuming that they are exposed to the background only, may therefore have 
produced spurious results. 

It is a pleasure to thank Jordi Miralda-Escude for enlightening discussions which greatly improved this work and 
Anthony Aguirre for a careful reading of the manuscript. I also thank the Institute for Computational Cosmology 
in Durham and the Institute of Astronomy in Cambridge, where part of this work was completed. Finally, I grate- 
fully acknowledge support from the W. M. Keck Foundation, from NSF grant PHY-0070928, and from Marie Curie 

^ The columns quoted are integrated over the systems (i.e., the columns of the individual Voigt profile components were summed) and 
the rates of incidence are for systems at redshift z = 2.5—3.5. We only used data with observed wavelength redwards of the quasar's hya 
emission line and with a redshift more than 4000 kms^^ bluewards of the redshift of the quasar. 
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